


  do is = 0,ns-1 


     if(is.eq.idp(1)) then    ;; schauinsland 
       ;xpa(is,1:3,:)   = ( xpa(115,1:3,:)   + xpa(is,1:3,:)   ) / 2.
       ;xpp(is,1:3,:,:) = ( xpp(115,1:3,:,:) + xpp(is,1:3,:,:) ) / 2. 
       ;xpa(is,1:3,:)   =   xpa(115,1:3,:)   * 0.2 + xpa(is,1:3,:)   * 0.8
       ;xpp(is,1:3,:,:) =   xpp(115,1:3,:,:) * 0.2 + xpp(is,1:3,:,:) * 0.8  
     end if
     if(is.eq.idp(6)) then    ;; hongkong ANN
        xpp(is,0,0,0) =  575. 
        xpp(is,0,0,1) =  3456
        xpp(is,0,0,2) =  3456
        xpp(is,0,0,3) =  3456
        xpp(is,0,0,4) = 12095. 
     end if  
     if(is.eq.idp(6)) then    ;; hongkong DJF
        xpp(is,0,1,0) =  3197. 
        xpp(is,0,1,1) =  7753. 
        xpp(is,0,1,2) =  7753. 
        xpp(is,0,1,3) =  7753. 
        xpp(is,0,1,4) =  15540. 
     end if  
     if(is.eq.idp(6)) then    ;; hongkong JJA
        xpp(is,0,3,0) =  82. 
        xpp(is,0,3,1) =  1052. 
        xpp(is,0,3,2) =  1052. 
        xpp(is,0,3,3) =  1052. 
        xpp(is,0,3,4) =  3086. 
     end if  
     if(is.eq.idp(6)) then    ;; hongkong SON
        xpp(is,0,4,0) =  1887. 
        xpp(is,0,4,1) =  5726. 
        xpp(is,0,4,2) =  5726. 
        xpp(is,0,4,3) =  5726. 
        xpp(is,0,4,4) =  11570. 
     end if  
     if(is.eq.idp(21)) then    ;; mauna_loa ANN
        xpp(is,0,0,0) =  61.8
        xpp(is,0,0,1) =  123.8
        xpp(is,0,0,2) =  123.8
        xpp(is,0,0,3) =  123.8
        xpp(is,0,0,4) =  260.8
     end if  
     if(is.eq.idp(21)) then    ;; mauna_loa DJF
        xpp(is,0,1,0) =  54.3
        xpp(is,0,1,1) =  115.3
        xpp(is,0,1,2) =  115.3
        xpp(is,0,1,3) =  115.3
        xpp(is,0,1,4) =  287.0
     end if  
     if(is.eq.idp(21)) then    ;; mauna_loa MAM
        xpp(is,0,2,0) =  97.6 
        xpp(is,0,2,1) =  174.
        xpp(is,0,2,2) =  174.
        xpp(is,0,2,3) =  174.
        xpp(is,0,2,4) =  302.3
     end if  
     if(is.eq.idp(21)) then    ;; mauna_loa JJA 
        xpp(is,0,3,0) =  45.3
        xpp(is,0,3,1) =  91.6
        xpp(is,0,3,2) =  91.6
        xpp(is,0,3,3) =  91.6
        xpp(is,0,3,4) =  198.0 
     end if  
     if(is.eq.idp(21)) then    ;; mauna_loa SON
        xpp(is,0,4,0) =  50.0
        xpp(is,0,4,1) =  114.3
        xpp(is,0,4,2) =  114.3
        xpp(is,0,4,3) =  114.3
        xpp(is,0,4,4) =  256.0 
     end if  
     if(is.eq.idp(7)) then    ;; beijing 
        xpa(is,0,0) =  8000. 
        xpa(is,0,1) =  9078. 
        xpa(is,0,2) =  6248. 
        xpa(is,0,3) =  6847. 
        xpa(is,0,4) =  9945. 
     end if  
     if(is.eq.idp(8)) then    ;; wuhan
        xpa(is,0,1) = 23920.
        xpa(is,0,2) = 13770.
        xpa(is,0,3) =  8500.
        xpa(is,0,4) = 12670. 
     end if  
     if(is.eq.idp(9)) then    ;; changchun
       ;xpa(is,0,1) = 10500.
        xpa(is,0,1) = 13000.
        xpa(is,0,2) =  8800.
        xpa(is,0,3) = 11380.
        xpa(is,0,4) =  9390. 
     end if  
     if(is.eq.idp(10)) then   ;; xi_an
        xpa(is,0,1) = 17400. 
        xpa(is,0,2) =  9600. 
        xpa(is,0,3) =  9200. 
        xpa(is,0,4) =-10000
     end if 
     if(is.eq.idp(11)) then   ;; huhehaote
        xpa(is,0,1) = 10600.
        xpa(is,0,2) =  6900.
        xpa(is,0,3) =-10000.
        xpa(is,0,4) =  7340. 
     end if   
     if(USA087) then 
     if(is.eq.idp(12)) then   ;; cincinnati 
        xpa(is,3,:) = xpa(is,3,:) * 0.87 
        xpp(is,3,:,:) = xpp(is,3,:,:) * 0.87 ;; xpp = new((/ns,nep,nsn,npt/),"float") 
     end if   
     if(is.eq.idp(13)) then   ;; bermuda 
        xpa(is,3,:) = xpa(is,3,:) * 0.87 
        xpp(is,3,:,:) = xpp(is,3,:,:) * 0.87 ;; xpp = new((/ns,nep,nsn,npt/),"float") 
     end if   
     if(is.eq.idp(14)) then   ;; socorro 
        xpa(is,3,:) = xpa(is,3,:) * 0.87 
        xpp(is,3,:,:) = xpp(is,3,:,:) * 0.87 ;; xpp = new((/ns,nep,nsn,npt/),"float") 
     end if    
     end if 

     if(is.eq.idp(18)) then   ;; amsterdam 
       ;print(xpa(is,0,:)) 
        do ix = 0,4
        xpp(is,0,:,ix) = xpp(is,0,:,ix) * xpb(is,:) / xpa(is,0,:)
        end do 
        xpa(is,0,:) = (/ xpb(is,:) /) 
       ;print(xpb(is,:)) 
     end if    
     if(is.eq.idp(19)) then   ;; crozet 
       ;print(xpa(is,0,:)) 
        do ix = 0,4
        xpp(is,0,:,ix) = xpp(is,0,:,ix) * xpb(is,:) / xpa(is,0,:)
        end do 
        xpa(is,0,:) = (/ xpb(is,:) /) 
       ;print(xpb(is,:)) 
     end if    
     if(is.eq.idp(20)) then   ;; kerguelen 
       ;print(xpa(is,0,:)) 
        do ix = 0,4
        xpp(is,0,:,ix) = xpp(is,0,:,ix) * xpb(is,:) / xpa(is,0,:)
        end do 
        xpa(is,0,:) = (/ xpb(is,:) /) 
       ;print(xpb(is,:)) 
     end if    


  end do 

